Effect of Resonant Neutrino Oscillation on TeV Neutrino Flavor Ratio 

from Choked GRBs 
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In the coUapsar scenario of the long duration Gamma-Ray Bursts (GRBs), multi-TeV neutrino 
emission is predicted as the jet makes its way through the stellar envelope. Such a neutrino signal 
is also expected for more general "failed" GRBs in which a putative jet is "choked" by a heavy 
envelope. If the Ve — >■ ^^^ neutrino oscillation parameters are in the atmospheric neutrino oscillation 
range, we show that the resonant oscillation of Ve <->■ v^^t can take place within the inner high 
density region of the choked jet progenitor with a heavy envelope, altering the v flavor ratio on its 
surface to "I>J_. : "I>J_ : $f,^ = 5 : 11 : 2. Considering vacuum oscillation of these neutrinos on their 
way to Earth, the final flavor ratio detected on Earth is further modified to either 1 : 1.095 : 1.095 
for the large mixing angle solution to the solar neutrino data, or 1:1.3:1.3 for maximal mixing among 
the muon and tau neutrinos in vacuum. 

PACS numbers: 98.70.Rz,96.40.Tv 



I. INTRODUCTION 

Long duration GRBs (LGRBs) are believed to be asso- 
ciated with deaths of massive stars [l[ . The evidence in 
support of such an origin includes associations of several 
LGRBs with Type Ic supernovae and the prevalence of 
star forming dwarf host galaxies associated with LGRBs 
[4I. Observationally, only a small fraction (< 10^"^) of 
core collapse SNe are associated with GRBs [3]. They 
correspond to those jets that break through the stellar 
envelope and reach a highly relativistic speed (Lorentz 
factor F > 100). Internal shocks are formed in the opti- 
cally thin regions, and gamma-rays are produced by syn- 
chrotron radiation and/or inverse- Compton scattering of 
Fermi accelerated electrons in these shocks. 

On the other hand, it is feasible to envision that a 
much larger fraction of core collapses may also launch a 
mildly relativistic jet from the central engine, but the jet 
never makes its way out from the envelope, due to either 
a smaller energy budget or a more extended, massive stel- 
lar envelope than in a GRB progenitor. In any case, both 
the successful and these "choked" jets can accelerate pro- 
tons to energy > 10^ GeV from the internal shocks well 
inside the stellar envelope. The interaction between these 
protons and the ^ 1 keV thermal X-ray photons emit- 
ted by the hot cocoon surrounding the jet would gener- 
ate multi-TeV neutrinos through photopion production 
[j|. For an individual GRB at rcdshift 2: ^ 1, the pre- 
dicted upward going muon event number is 0.1-10 in a 
km'^ detector [J]. There can also be neutrino production 
due to pp and pn collisions involving relativistic protons 
from the buried jet and the thermal nucleons from the 
jet and the surrounding, which can produce more abun- 
dant neutrinos for the presupernova stars with a heavy 
envelope ^ . The detection of low luminosity (LL) GRBs 
(such as GRB 980425 and GRB 060218 0) suggests that 
the event rate of gamma-ray dim core collapses is much 
higher than those associated with high luminosity GRBs 



Q. It is conceivable that the gamma-ray "dark" choked 
GRBs are even more abundant, and would contribute 
more to the high energy neutrino background. 

Several other neutrino mechanisms have been dis- 
cussed in the literature. The internal shocks that power 
the prompt gam ma-ray emission can produce ^100 TeV 
neutrinos jsMlOJ , which should be lagged behind the TeV 
neutrinos. Inelastic collisions between decoupled protons 
and neutrons during the acceleration of the fireball can 
power a multi-GeV neutrino signal [ll| , but the predicted 
flux level is below the atmospheric neutrino background, 
and hence, difficult to detect. The neutrinos produced 
in the early afterglow phase have a very high energy (~ 
EeV) [141 and are not optimized for detection with the 
present day neutrino detectors. The multi-TeV neutrino 
signal discussed here could be above the atmospheric 
background, and may be detectable with km'^ detectors 
for some nearby sources. 

The main source of high energy neutrinos is the decay 
of charged pions, which leads to the neutrino flux ratio 
at the production site $°^ : $°^ : $°^ = 1 : 2 : ($°^ 
corresponds to the sum of neutrino and anti-neutrino flux 
for the flavor a at the source). The vacuum oscillation 
of these neutrinos on their way to Earth would make the 
observed ratio to 1 : 1 : 1. This applies to low energy 
neutrinos including the TeV neutrinos discussed in this 
paper. For high energy neutrinos above ~1 PeV, muon 
energy is degraded before decaying to low energy neutri- 
nos so that high energy neutrinos will be absent. The 
neutrino flux ratio at the source is modified to : 1 : 0, 
which is further modified to 1 : 1.8 : 1.8 at Earth after 
vacuum oscillation is taken into account |13| . 

Another possibility of modifying the neutrino flavor 
ratio is the resonant conversion of neutrinos from one 
flavor to another due to the medium effect. Such an effect 
is known to be important for solar neutrinos [14J, and 
has been discussed for hot plasma in the early universe 
[l5| . supernova medium [ig, the GRB fireball uE 



and jet [2i|. Here we show that for choked GRBs, the 
multi-TeV neutrinos discussed by Meszaros & Waxman 
Q could undergo resonance oscihation in the high density 
core (typicaUy He core) of the presupernova star, if the 
neutrino oscillation parameters are in the atmospheric 
neutrino oscillation range. This would alter the neutrino 
flavor ratio escaping from the stellar envelope, and hence, 
the eventual detected flavor ratio from Earth. 



II. NEUTRINO OSCILLATION IN THE 
STELLAR ENVELOPE 

As a mildly relativistic jet makes its way through the 
stellar envelope, internal shocks can develop and can ac- 
celerate protons to energy ~ 10^ GeV. These protons 
would interact with the '^ keV thermal X-ray photons 
to produce ~ 5 TeV neutrinos via the process p + 7 ^■ 
A+ -^ n + TT'^ ^ n + fj,~^ + Vf^ ^ n + e^ + Vfj, + i^e + i^fj.- 

Depending on the initial mass and metalicity, the pre- 
supernova star can have different compositions with dif- 
ferent radii. The LGRB progenitors (Type Ic SN) have 
lost the H and most of the He envelopes before explo- 
sion. They are too small to have an interesting neutrino 
oscillation signature. The choked jet progenitors, on the 
other hand, can retain the He envelope (Type lb SNe) 
and even the H envelope (Type H SNe) . These presuper- 
nova stars are favorable for TeV neutrino production and 
neutrino oscillation. It is believed that the putative jet 
is launched along the rotation axis where the centrifugal 
support is the least, and is powered by either vi> annihila- 
tion or through some electromagnetic processes. Without 
exploring the details of jet dynamics, here we treat the 
jet parametrically. We assume that the jet is developed 
and the TeV neutrinos are produced at a radius rj <C -R* , 
where i?* is the radius of the star. 

Depending on the energy of the propagating neutrino 
and the nature of the background, neutrinos can interact 
with the background particles via charge current (CC) 
and neutral current (NC) interactions. For a neutrino 
energy below E^, ~ M^/2me — 10^ GeV, an electron 
neutrino can have both GG and NG interactions with 
the normal matter, whereas muon and tau neutrinos can 
only have NC interactions. The effective potential of NG 
interactions is the same for all active neutrinos. Since 
oscillation depends only on the potential difference, for 
active-active oscillations the NG contributions cancel out. 
So only the GG contribution to the neutrino potential, 
V — y^GpNe, is responsible for neutrino oscillation in 
the medium, where Gp is the Fermi coupling constant 
and Ne is the electron number density in the medium. 
For anti- neutrinos N^ is replaced by — iVg. Thus for the 
process Ve ^ ^ai,t, the neutrino potential is ^/2GFNe, 
while for the process !/^ O Vr it vanishes. 

Here we consider the simplified picture of two mixed 
flavor states v^ and 1/^(1/.^) with the vacuum mixing an- 
gle and mass square difference Am?. In an uniform 
medium the evolution of the flavor states are governed 
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where A — Am? /2Eu, V is the potential difference be- 
tween Vu^ and Vu^, (i.e. V = V^^ — Vu^) and Eu is the 
neutrino energy. The transition probability as a function 
of distance £ is given by 
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Once the neutrinos are produced due to pion decay at 
a point Tj ^ i?^, they will propagate away from the 
star where the medium effect can be substantial. If 
the density of the medium is such that the condition 
^/2GFNe = A cos 20 is satisfied, then resonant conver- 
sion of neutrinos from one flavor to another with maxi- 
mum amplitude can occur. For anti-neutrinos the reso- 
nance condition can never be satisfied (for normal neu- 
trino mass hierarchy). Although the oscillation process 
De -f-)- Vfi^T can take place, it will be suppressed. 

The critical density for resonance is called the reso- 
nance density. For 5 TeV neutrinos, it reads 
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where we have Am? in units of eV^ and E^i2,y in units 
of 10^^'^ eV. The resonance length is 
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Define the stellar radius rpj as the radius at which the 
local density is pR. The first condition for resonant os- 
cillation is ^R, < tr. 

If the resonance region is wide enough the transi- 
tion can be total. We can define a resonance width for 
which the amplitude of the probability can be 1/2 in- 
stead of unity. In this case the width can be given as 
r = 2Am^ sin 20. This corresponds to a length scale 



Str = 



2 tan 20 

I ATe dr I 1 
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For (5rR > £r, there can be enough time for i/g to stay in 
the resonance region to convert into v^ivT-). This is the 
second condition for significant resonant oscillation. 

In order to evaluate both conditions, one needs to know 
the matter density profile in the stars (which determines 
tr and dNe/dr). The density profile of a presupernova 
star is difficult to probe observationally. Numerical mod- 
els predict a decreasing density with radius. If convective 



mixing is not important, there is a sharp decrease in den- 
sity beyond the He core with radius thc ~ 10^^ cm and 
the local density pHc ^ 10^^ g cm^'^. If convective mix- 
ing is important, there is no abrupt transition, and the 
density profile may be roughly described in the analytical 
formg Q 



r 



1 



(7) 



The parameters i?* and po depend on the type of the 
star. For example, a blue supergiant (BSG) model for SN 
1987A gives i?* = 3 x 10^^ cm and po = 3 x lO"'"^ g cm^ 



24j . In some models, the He core can extend to 10^^ cm. 



and the H envelope can extend to 10^'^ cm [J, l5| 

The exponent n = 3, 3/2 correspond to the radiative 
and convective envelope, respectively. In general it can 
vary between 2 and 3 for different numerical models. The 
condition Sr^ > £yi can be re-written as a requirement 
to n (where r ~ €r has been adopted which is relevant 
for resonant oscillation) 



n < 2tan26' 1 






(8) 



The resonance length l-^ depends on neutrino oscillation 
parameters and neutrino energy. Apparently, if ^r, > 
i?^, the requirement to n (n < 0) is unphysical, and no 
neutrino oscillation is expected inside these stars. On the 
other hand, if ^r <C i?*, the constraint n < 2 tan 261 may 
be satisfied in some stars for some oscillation parameters. 
In order to evaluate whether the neutrino oscillation 
conditions are satisfied, one needs to know the neutrino 
oscillation parameters in matter. Experimentally these 
are inaccessible. Only oscillation parameters from the so- 
lar and atmospheric neutrino experiments are available. 
For small neutrino mixing angles, the mixing matrix is al- 
most diagonal and each flavor eigenstate nearly overlaps 
with one of the mass eigenstate. One may then associate 
Ve to vi, Vfi to 1^2, and Vr to V3. While the solar neu- 
trino oscillation parameters are relevant to Vg — >■ u^ os- 
cillations, the atmospheric neutrino oscillation data cor- 
respond mostly to j/^ — > i^r oscillations, and the corre- 
sponding neutrino parameters are 023 and Am^^ , respec- 
tively. It is also very possible that the physical prop- 
erties of neutrinos in a medium could be different from 
their vacuum values. For example, a neutrino can acquire 
mass due to its interaction with the background particles 
even if we consider it massless in vacuum. Similarly the 
mixing properties in matter may not follow the vacuum 
pattern as measured. Nonetheless, since the oscillation 
parameters of the solar and atmospheric experiments are 
the best measured, in the following we test whether these 
parameters may allow neutrino oscillation to happen in 
the progenitor stars of choked GRBs. We do not take 
it for granted that any of these parameters are operat- 
ing in the oscillation process i^e ^ J^/^.t in the choked 
fireball, but just take the only experimentally available 
parameters to test the conditions for prominent oscilla- 
tions. Similar analyses have been carried out before to 



evaluate the possible oscillation effect in GRB fireballs 

The SNO salt phase solar neutrino data, combined 
with the KamLand reactor antineutrino results, constrain 
the neutrino oscillation parameters in the regime 6 x 
10-5 eV^ < Am2 < IQ-^eV^ and 0.64 < sin^ 26* < 0.96 
[25^, with the best fit parameters Am^ ~ 7.1 x lO'^ eV^ 
and sin^ 29 - 0.69 with 99% confidence level. The best 
fit values give pr,sno — 5.2 x lO^^g cm^^^iJ^Jj.T and 
%,SNO — 2.1 X 10^^cmi?i,,i2.7- We can see that £u is 
larger than i?* of a typical BSG, suggesting that resonant 
oscillation would not occur for these neutrino oscillation 
parameters. 

On the other hand, the atmospheric neutrino oscilla- 
tion parameters reported by the SK Collaboration are in 
the range 1.9 x lO-^eV^ < Am^ < 3.0 x lO'^eV^ and 
0.9 < sin^ 26 < 1.0 at a 90% confidence level ^, which 
corresponds to the oscillations of mostly muon neutri- 
nos to tau neutrinos. If we assume that these param- 
eters apply to j/g — > h'^.j- oscillations in matter, we can 
get the following constraint. We consider the good fit 
point Ato^ -- 2.5 x 10"^ eV^ and sin^ 26^ ~ 0.9, and 
get PR,SK :i; 1.0 X 10~3 gcm-3 E~l^^ and £r,sk ^ 
5.2 X 10^^ cm Ey^i2j. For the nominal BSG model 
discussed in this paper, the stellar radius at which the 
density is pr.sk is ?'r — 7.1 x 10^^ cm for n — "i and 
Eu,i2.7 = 1- We can see that the condition ^r < tr 
is satisfied for the typical neutrino energy E,y = 5 TeV. 
By using the value of tan 29 from the SK neutrino data 
and £ji = %,sk, we obtain from Eq. ([8]) n < 4.96. Known 
stellar models have n between 2 and 3. This suggests that 
the second condition is also fully satisfied. We conclude 
that resonant oscillation of multi-TeV neutrinos can oc- 
cur within a nominal BSG progenitor for the neutrino 
oscillation parameters inferred by the atmospheric neu- 
trino data. A similar analysis suggests that the same 
conclusion applies to other BSG progenitors or He stars 
with extended envelopes (with R^, up to 10"'^^ cm), but 
does not apply to typical He stars (with R^, — 10^^ cm), 
or other more compact stars. Since GRB observations 
favor associations with Type Ic SNe (for which the He 
envelope is mostly stripped off), GRBs are not preferred 
sources for TeV neutrino oscillation. Instead, we identify 
choked GRBs, especially those with a heavy envelope, as 
interesting sources for resonant neutrino oscillation. 

For a full oscillation, v^ can oscillate to i'^ and to ^'r 
with equal probability but i/^ can oscillate only to i^e- 
On average we can have 1/3 survival probability for i^e, 
i/^ and Vr for each v^ oscillation, but have 1/2 survival 
probability for Vg, v^ for the v^ -f-J- v^ resonant oscilla- 
tion. The p7 process also produce D^ which does not 
resonantly oscillate. Putting together, on the surface of 
the presupernova star the survival probability of each fla- 
vor (both neutrino and anti-neutrino) is in the ratio of 
(l + 5) • (i + 5 + 1) ■ 5 = 5 : 11 : 2. Above rR (mostly in 
the H envelope), the density is much lower than the res- 
onance density, and no further reconversion of neutrinos 
can take place. So the ratio 5 : 11 : 2 is final neutrino 



flavor ratio escaping from the star. This ratio is notably 
different from the nominal 1:2:0 ratio for multi-TeV 
neutrinos without considering the resonant oscillation ef- 
fect. 



III. NEUTRINO OSCILLATION IN VACUUM 

Since these choked GRB neutrino sources are typically 
at large distances, the TeV neutrinos escaping from the 
star would undergo vacuum oscillation on their way to 
Earth. The neutrino flux for a particular flavor a on 
Earth is given by 



*. 



E^" 



$ 



I'/j' 



(9) 



where <&* signifies the flux of i/^ at the surface of the He 
envelope after resonant oscillation, and Pap corresponds 
to the oscillation probability from h'a to i^p in vacuum. 
For the matter effect on resonant oscillation in the stellar 
envelope, we have applied a two-flavor neutrino analysis. 
This is because no Vr 's are generated in the P7 process on 
one hand, and the two-flavor neutrino analysis is simpler 
on the other hand, as it depends only on one mass square 
difference. The limited size of the oscillation baseline (the 
stellar envelope) also makes the three-flavor oscillation 
effect unimportant. Such a two-flavor neutrino analysis 
has been applied in most previous resonant oscillation 
discussions [l8l[22|. 

When discussing the vacuum oscillation effect along 
a long base line from the source to Earth, one needs 
however fully take into account the three-flavor neutrino 
oscillation effect. This is demanded by the combined 
analyses of both the solar and the atmospheric neutrino 
anomalies. For the best fit to the SNO data from the 
large mixing angle (LMA) solution one can take the mix- 
ing angles 0i2 = 34° ± 2.5°^23_= 45° ± 6°, fi^ia = 0° ± 8° 
and the Dirac phase 6 = [ij,l23]- This gives Pee — 0.57, 



efi 



Pe 



0.215 and P„„ = P^r = Prr ^ 0.393. In- 



serting these probabilities in Eq. © and considering the 
error of 0i2, the flux ratio at Earth is 1 : (1.095 ±0.012) : 
(1.095 ± 0.012). On the other hand, if we consider the 
maximal mixing among the i^^ and Vt in vacuum, then 
the Ve oscillation to t^^ and i^r is largely suppressed. One 



1, ^MM = Prr = P,.r = 1/2, and aU 



can then have Pe, 
other transition probabilities are negligible. Using these 
oscillation probabilities, we obtain the flux ratio at Earth 
as 1:1.3:1.3. 



IV. DISCUSSION 

GRBs have a wide redshift distribution (from z = 
0.0085 to z = 8.2). Observations suggest that the nearby 
low-luminosity GRBs have a local event rate^ 100 time 
higher than that of high-luminosity GRBs [7]. If GRB 
jets become progressively successful in progressively rarer 



progenitors, it is conceivable that there could be even 
more choked GRBs in the nearby universe. Although the 
gamma-ray luminosity becomes progressively smaller as 
the envelope becomes progressively heavier, the TeV neu- 
trino luminosity may not decrease, and could even follow 
an opposite trend. Assuming that the choked GRB pro- 
genitor has a local event rate similar to that of LL-GRBs, 
i.e. TZ ~ 200 Gpc"'^ yr~^, one would expect ~ 14 neu- 
trino bursts (without gamma-ray counterpart) per year 
all sky at z < 0.1. If the TeV neutrino luminosity of these 
events are similar to those of successful GRBs [4, 5], then 
each event would have 100s of TeV neutrinos detected by 
a km-^ detector such as IceCube. Such a possibility is al- 
ready ruled out by the current upper limits placed by the 
IceCube observations. This suggests that either there are 
not that many nearby choked GRBs, or that the choked 
GRBs are not as neutrino- luminous as predicted [4] . Go- 
ing to the conservative extreme, i.e. if the neutrino lu- 
minosity is correlated with gamma-ray luminosity, then 
the detected event rate for these nearby neutrino burst 
sources would be 0.001-0.01 neutrinos per event in a km^ 
detector. This is essentially impossible for IceCube to de- 
tect individual sources. The real detected neutrino event 
rate may be between these two extreme values. IceCube 
or a similar detector would be able to detect these neu- 
trino bursts or to place even more stringent upper limits 
in the near future. 

If the nearby neutrino bursts are bright enough, the 
deviation of the observed neutrino flavor ratio from 1:1:1 
may be tested by IceCube or similar detectors. The fla- 
vor ratios can be in principle deduced from the relative 
rates of showers, muon tracks, and the unique tau lep- 
ton induced signals [28]. The possibility of detecting a 
tau signal by IceCube is low, especially in the multi- 
TeV energy range. On the other hand, IceCube can 
distinguish between shower-like events and the /.t-track 
events, although it is hard to identify i^e and Vr through 
their electromagnetic and hadronic showers. Nonethe- 
less, assuming a flavor-independent neutrino spectrum 
and i/^ — v-r symmetry (as is the case of our two pre- 
dicted ratios), Ve fraction may be extracted from the 
measured Muon to shower ratio [2g. The 10% differ- 
ence in the flavor ratio reduces the Ve fraction from 
1/3 to 0.313 (for flavor ratio 1:1.095:1.095). This cor- 
responds to a slight increase of Muon to shower ratio. 
With the uncertainty (20%) for the nominal diffuse flux 
{El^dN^JdE^^^ = 10"^ GeVcm-2s-i for one year) 
adopted in [28j] , the small change in Muon to shower ratio 
may not be differentiated. If nearby neutrino bursts are 
bright enough, the flux would be increased and the un- 
certainty reduced significantly. This would make a bet- 
ter case for detecting the flavor ratio change. For the 
1:1.3:1.3 ratio, the Vg fraction is reduced from 1/3 to 0.28, 
making the Muon to shower ratio as high as ~ 3.5 (as 
compared to ~ 3 for 1:1:1). The effect may be detectable 
for the putative bright neutrino burst events discussed 
above, even if they may be very rare. 

Since the parameters (£r, pr, tr, and 5r^) all de- 



pend on E^, we expect the flavor ratio would also de- 
pend on neutrino energy. This aspect has been extensive 
discussed in [23|- 
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